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We report a detailed study of the pairing resonance via tunneling density of states in ultra-thin 
superconducting Al films in supercritical magnetic fields. Particular emphasis is placed on effects of 
the perpendicular component of the magnetic field on the resonance energy and magnitude. Though 
the resonance is broadened and attenuated by H± as expected, its energy is shifted upward linearly 
with H±. Extension of the original theory of the resonance to include strong perpendicular fields 
shows that at sufficiently large H± the overlap of the broadened resonance tail with the underlying 
degenerate Fermi sea alters the spectral distribution of the resonance via the exclusion principle. 
This leads to the shift of the the resonance feature to higher energy. 

PACS numbers: 74.25.Qt,74.25.Op,74.40.+k 



Though superconductivity is one of the most ubiquitous and extensively studied phenomenon in condensed matter 
physics, little is known about how an attractive electron-electron interaction mode, whether mediated by phonons or 
some more exotic coupling, evolves from a non-coherent fluctuation channel to a full blown superconducting instability. 
Because the formation of stable Cooper pairs usually coincides with a precipitous collapse into a macroscopic quantum 
condensate, it is difficult to extract properties of individual Cooper pairs in the superconducting phase. Likewise there 
are few systems in which a Cooper pair precursor can even be identified in the normal state Q,y>Q, much less studied. 
Here we show that spin-paramagnetically limited superconducting Al films provide an unique opportunity to directly 
probe the energetics and dynamics of the Cooper pair precursor near the threshold of a well-defined superconducting 
instability associated with the first-order spin-paramagnetic (SP) transition. Using electron tunneling spectroscopy, 
we demonstrate that the precursor, which is manifest as a finite bias anomaly in the tunneling spectrum, exhibits an 
unexpected energy shift associated with the orbital response to the perpendicular component of the applied magnetic 
field. 

Although the spin-paramagnetic transition was conjectured over forty years ago y| la , the dynamics and complete 
phase diagram of the transition have only been mapped out in the last decade |ti ItI W I'l l 1(1 1 1 1| . This is, in part, due 
to the fact that the first-order nature of the transition and the corresponding tricritical point (Tth ~ 600 mK in Al) 
are extremely sensitive to and quickly extinguished by spin-orbit (SO) scattering. Consequently, many of the salient 
features of the SP transition such as hysteresis 0, state memory 0], avalanches [l^, and the finite bias tunneling 
anomaly require a true spin-singlet ground state. This precludes the observation of such effects in high atomic 
number superconductors in which there is a strong spin-orbit scattering rate, such as there is in Pb and Sn, for 
instance [lj,ll5j. Indeed, there are only two well documented spin-singlet superconductors, Be and Al. 

In practice, the SP transition can be probed by applying a magnetic field parallel to the surface of a film that is 
sufficiently thin (< 5 nm) so as to suppress the Meissner response. The field thus fully penetrates the film, and the 
superconducting phase remains unperturbed until the Zeeman splitting becomes of the order of the superconducting 
gap. At that point a first order phase transition to the paramagnetic normal state occurs 0,0- The spin symmetry of 
ultra-light BCS superconductors is not only manifested in the SP transition itself but in the paramagnetic normal state 
as well. In particular, an unexpected pairing resonance (PR), clearly associated with virtual Cooper pair formation, 
was recently discovered in the paramagnetic normal state density of states (DOS) spectrum of Al and Be films 0,^1- 
The resonance can be observed in high field tunneling by virtue of the fact that it is a spin-singlet mode riding on top 
of a paramagnetic background. In general terms its energy is the sum of the Zeeman energy of the anti-aligned electron 
in the virtual Cooper pair minus the effective binding energy gained by forming the pair. A rigorous nonperturbative 
analysis by Altshuler et. al 0, showed that in parallel field the resonance occurs at an energy that is universal 



2 



for OD(grain), lD(wire), and 2D(film) systems, 

E+ = (E z + fi)/2, (1) 

where Q = — A 2 , E z — ql^bH is the Zeeman energy, H the magnetic field in Tesla, g^ the Lande g-factor, 

Hb the Bohr magneton, and A G the zero temperature, zero field, superconducting energy gap. Though Eq. Q was 
found to be in reasonably good agreement with parallel field experiments |l9l |. a comprehensive comparison between 
experiment and theory was never made. In particular, though it is known that the resonance becomes more stable 
with decreasing film conductance, its evolution in the limit that the dimcnsionlcss conductance g — h/e 2 R — > 1, where 
R is the sheet resistance of the film, remains an open question. Furthermore, the theory was developed under the 
assumption that the Zeeman energy dominates the spectrum. This assumption will be violated in a tilted field if the 
Cooperon cyclotron energy MIh = ADeH± ~ E z , where D is the electron diffusivity and H± is the perpendicular 
component of the field. In this Letter, we take advantage of improved sample quality to investigate the PR in low g 
films in the presence of large perpendicular magnetic fields. We show that contrary to expectation, the PR energy in 
tilted field is not solely a function of the Zeeman splitting, but, in fact, has a contribution that is proportional to the 
perpendicular component of the field. 

The films were grown by e-beam deposition of 99.999% Al onto fire polished glass microscope slides held at 84 
K. The depositions were made at a rate of ~ 0.1 nm/s in a typical vacuum P < 3 x 10~ 7 Torr. A series of films 
with thicknesses ranging from 2 to 2.9 nm had a dimcnsionlcss normal state conductance that ranged from g = 5.6 
to 213 at 100 mK. After deposition, the films were exposed to the atmosphere for 0.5-4 hours in order to allow a 
thin native oxide layer to form. Then a 9-nm thick Al counterelectrode was deposited onto the film with the oxide 
serving as the tunneling barrier. The counterelectrode had a parallel critical field of —2.7 T due to its relatively large 
thickness, which is to be compared with H c n ~ 6 T for the films. The junction area was about 1 mmxl mm, while 
the junction resistance ranged from 10-100 kfi depending on exposure time and other factors. Only junctions with 
resistances much greater than that of the films were used. Measurements of resistance and tunneling were carried out 
on an Oxford dilution refrigerator using a standard ac four-probe technique. Magnetic fields of up to 9 T were applied 
using a superconducting solenoid. A mechanical rotator was employed to orient the sample in situ with a precision 
of - 0.1°. 

Two typical 60-mK tunneling spectra are shown in Fig.^in which the parallel field is above and below H c n . All of the 
tunneling data presented in this paper were made with the counter-electrode in the normal state. In this configuration 
the low temperature tunnel conductance is proportional to the DOS of the film [2jj. In the superconducting state, 
the BCS DOS is Zeeman-split by the applied magnetic field as can be seen in the 4.5-T curve of Fig. CD As one further 
increases the field, a first-order transition to the normal state occurs near H c \\ = A Q / \fgZp,B |3. M. l2lj. The normal 
state spectrum, represented by the 5.7 T curve in Fig. ^ displays a logarithmic depletion of states near the Fermi 
energy, commonly known as the zero bias anomaly (ZBA). It is now well established that the logarithmic ZBA is 
associated with electron-electron interactions in two-dimensional disordered systems j^. The finite bias anomalies 
in the 5.7-T curve are due to the PR and represent a window into the processes that ultimately lead to the formation 
of stable Cooper pairs at lower fields. 

In Fig. we plot the resonance feature in a g = 5.6 sample at a variety of tilt angles 9, where 8 is the angle 
between the field and the plane of the film. The field was held constant in this data set so that the perpendicular field 
component was H± = H sin 6*. Two features of the data should be noted. First, though the PR is attenuated and 
broadened by H±, it survives to much larger tilt angles in the lowest conductance films. This is clear in the inset of 
Fig. [21 where the relative depth of the PR is plotted as a function of H± for films of different g. Second, the resonance 
moves to higher energy with increasing tilt angle. 

The measurement of PR in tilted magnetic fields allows one not only to test the effect of perpendicular field, but 
also extend the measurements to fields well below the spin-paramagnetic critical field, H c it. Upon rotation of a few 
degrees out of parallel the first-order spin-paramagnetic transition becomes second-order and the critical field rapidly 
begins to decrease until finally reaching H C 2 at 8 — 90°. This latter behavior can be seen in the inset of Fig. EI In the 
main panel of Fig. |3 we plot the resonance voltage as a function of applied field for the g = 5.6 sample. The solid 
symbols represent parallel field data which, of course, must terminate at H c \\ due to the fact that the film undergoes 
a first-order transition to the superconducting phase at this field. The solid line represents a fit to the parallel field 
data using Eq. Q), where only the g-factor was varied, with a best fit value gh = 1-61 0. The gap value used in 
Eq. was taken from measurements in the superconducting phase, A G /e = 0.43 mV. The open symbols in Fig. [3] 
represent measurements made just above the critical field at finite tilt angles up to 8 — 90°. Note that the tilted field 
data are linear in H well below if c || and that this behavior is inconsistent with Eq. 

The open symbol data points in Fig. [3] were obtained by changing both the field magnitude and the tilt angle. A 
more straightforward characterization of the field dependence of the resonance can be obtained by varying the field 



3 



at a preset tilt angle. Shown in Fig. 0] is the field dependence of the resonance energy at tilt angles 6 = 0°, 23°, 35°, 
and 90° . Not only is the linearity of the parallel field dependence preserved in tilted field, but the slope of the field 
dependence grows with increasing tilt angle. 

We have extended the original analysis of the resonance to include the case where the perpendicular magnetic field 
is such that the Cooperon cyclotron energy is much larger than the characteristic width of the resonance 

A 2 

Ml H = AeDHsinO > W 2 = —^k- (2) 

Agil 

The above limit violates an approximation used in the derivation of Eq. JJ| that requires the Cooperon momentum 
transfer to be small. In the parallel field case, the transition from the superconducting phase to the paramagnetic 
phase occurs at the Clogston limit E z ~ y/g~[Ao', consequently the evanescent pair energy fl lies well above the Fermi 
energy, fl > A„. Assuming that the momentum transfer is small, i.e. W2 <§C A Q , one can neglect the presence of 
the degenerate Fermi sea below the pair, since E + > CI. In the presence of a perpendicular magnetic field, however, 
Landau quantization of the electron motion leads to a momentum transfer (in energy units) of the order of hfln, and 
when this becomes sufficiently large the exclusion principle can no longer be ignored. When taking into account the 
Fermi sea, the position of the resonance feature is shifted to higher energies, as we explain by the diagram in the 
upper inset of Fig. 0] On the left we show qualitatively the anomaly in the density of states for Kin <C W% ■ the 
anomaly is centered at E + and its width is given by W2; since W% <C E + , the presence of the Fermi sea (gray area) 
is not felt by the excitations determining the anomaly. On the right, the dashed line shows the anomaly as it would 
be without the Fermi sea for MIr 3> W2', a finite fraction of excitations contributing to the anomaly would now be 
located below the Fermi energy. The solid line represents the "true" shape of the anomaly. As the exclusion principle 
suppresses the contribution of the excitations below the Fermi energy, the profile of the anomaly becomes asymmetric 
and the observed minimum is shifted to higher energies; moreover we expect that the larger Qh is, the larger the shift 
becomes, in agreement with the quantitative result given below. 

As will be described in detail elsewhere [2^] , the total correction to the density of states can be written as the sum 
of two terms, 8v = 8i>\ + 8V2, where the first is the perpendicular field result previously obtained in Ref. ^3 an d the 
second is the new contribution of the present analysis. We reemphasize that the minimum of 8v\ is always centered at 
E+, for any tilting angle. While the full expression for 8v2 is somewhat lengthy and beyond the scope of this Letter - 
it is obtained by a perturbative (one- loop) calculation, justified by the condition in Eq. (J2J), which includes subleading 
contributions - it is this term that shifts the apparent position of the resonance to higher energy. By linearizing the 
energy derivative of 8v, we obtain the following approximate formula for the position V* of the resonance in tilted 
field: 

V* = E + /e + a4DH$ai6, a = (~) - 0-055 (3) 

where ip is the digamma function. This generalizes the result of Rcf. and reduces to it in the limit H± — 0. The 
E + term in Eq. |J3J is well approximated by a linear field dependence, see the dashed line Fig. [31 and has no angle 
dependence. Consequently, the contribution of the orbital term can be extracted from the linear fits to the tilted 
field data in Fig. 21 (solid lines). The slopes values obtained from these fits are plotted as a function of sin (9 in the 
lower inset of Fig. 0| The dashed line in the inset has a slope of 0.38fi,/m which by Eq. should be equal to aAD. 
This gives a corresponding diffusivity D ~ 1.7h/m which is in excellent agreement with the value of the diffusivity 
obtained from the film conductance D = g/Airhft ~ lAh/m, where v is the renormalized density of states and t ~ 1 
nm is the film thickness after accounting for surface oxidation. 

In conclusion, we have studied the effect of perpendicular field on pairing resonance in the low temperature para- 
magnetic normal state of ultra-thin Al films. Tilting the magnetic field away from parallel orientation broadens the 
resonance thereby producing an overlap with the underlying Fermi sea. The exclusion principle effectively produces 
a "boundary condition" that truncates the low side of the resonance tail, resulting in an apparent upward shift in the 
resonance energy. Since the resonance represents a microscopic window into the fundamental quantum processes by 
which a stable Cooper pair is formed out of a specific interaction channel, a more comprehensive understanding of 
its dynamics has far reaching implications. In particular, further studies in granular Al films with significantly lower 
conductance may reveal a dimensionality crossover from 2D to 0D as g — > 1, as well as illuminate the role of e-e 
interaction effects in the evolution of the resonance. 

We gratefully acknowledge enlightening discussions with Igor Aleiner, Ilya Vekhter, Dana Browne, and Nandini 
Trivedi. This work was supported by the National Science Foundation under Grant DMR 02-04871. 
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FIG. 1: Tunneling conductance of a g — 12 film at 60 mK in superconducting and normal states, where parallel fields are 4.5 
T (a) and 5.7 T (b), respectively. In the superconducting state, the Zeeman split of BCS DOS manifests itself as two peaks of 
tunneling conductance on either side of V = 0. 
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FIG. 2: Pairing resonance at several field orientations for a g = 5.6 film at 60 mK in a field of 5.9 T. 6 = corresponds to 
field parallel to the film plane. The vertical dashed line marks the position of the parallel field resonance. Inset: the relative 
depth of the PR as a function of the perpendicular component of magnetic field for several film conductances: Circle g = 58, 
Square g — 26, Diamond g — 12, Triangle g — 5.6. 
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FIG. 3: The position of the PR as a function of magnetic field for the g = 5.6 film at 60 mK. The vertical ar- 
row marks the parallel critical field. Solid symbols correspond to parallel field. The open symbols were taken at 9 = 
8°, 12°, 16°, 19°, 23°, 29°, 45°, and 90° (from right to left). The dashed line is provided as a guide to the eye. The solid 
line is a fit of Eq. Q to the parallel field data. The inset shows the critical field as a function of tilt angle. 
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FIG. 4: The position of the PR as a function of magnetic field for the g — 5.6 film at 60 mK at several different field 
orientations: circle 9 — 0° , up triangle 9 — 23° , diamond 9 = 35° , and down triangle 9 = 90° . The solid lines are linear fits to 
the data. For comparison, the Zeeman voltage V z = E z /e is plotted as a dashed line, taking §l = 2. Upper inset: Width and 
position of the resonance for small (left) and large (right) perpendicular field as explained in the text. Lower inset: the fitted 
slope dV* /dH in units of h/m from the data in the main panel plotted as a function of sin^. The dashed line is a linear fit to 
the data and has slope 0.38. 



